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Abstract. In this paper, we are interested in tlie nonlinear Schrodinger problem —Aii + 
Vu = \u\f^^u submitted to the Dirichlet boundary conditions. We consider p > 2 and 
we are working with an open bounded domain Q. C M*^ (N > 2). Potential V satisfies 
max(l',0) e L'^I^(Q.) and min(y,0) e L+°°(i2). Moreover, -A + V is positive definite 
and has one and only one principal eigenvalue. When p ~ 2, we prove the uniqueness 
of the solution once we fix the projection on an eigenspace of —A + V. It implies partial 
symmetries (or symmetry breaking) for ground state and least energy nodal solutions. In 
the litterature, the case V = has already been studied. Here, we generalize the technique at 
our case by pointing out and explaining differences. To finish, as illustration, we implement 
the (modified) mountain pass algorithm to work with V negative, piecewise constant or not 
bounded. It permits us to exhibit direct examples where the solutions break down the 
symmetries of V . 



1. Introduction 

Let N>2, p>2, X>0 and an open bounded domain Q. C . We study the nonlinear 
Schrodinger problem 

-Au{x)+V{x)u{x)^X\u{x)\P'^u{x) (1) 

submitted to the Dirichlet boundary conditions (DBC). We are interested in the symmetry 
of solutions. 

When V belongs to L^/^(i2), the solutions can be defined as the critical points of the 
energy functional 



2 Ja p Ja 

Clearly, is solution. Concerning other solutions, if we assume that —A + V is posi 
tive definite and V := min(y,0) G L+°°(n), then the norm ||m|P = /fj|VMp + yM^ defined 
on Hq (n) is equivalent to the traditional norm ||m||?^i — /s^IVm^ (see Proposition 2.1 1. By 



working in the same way as in [8], it directly implies the existence of ground state solutions 
(g.s.) and least energy nodal solutions (l.e.n.s.); i.e. one-signed (resp. sign-changing) solu- 
tions with minimal energy. These solutions are characterized as minima of Sp respectively 
on the (resp. nodal) Nehari set 



Me//o'(n)\{o} I / \Vu\^+vi? = x f \u\p 

Jn Jo. 



(resp. :— {m : e ^p})- The Morse index is 1 (resp. 2). 

In this paper, we study the structure of these two types of solutions. We verify whether 
they are odd or even with respect to the hyperplanes leaving V invariant (i.e. V respects 



1991 Mathematics Subject Classification, primary 35J20; secondary 35A30. 

Key words and phrases. Nonlinear Schrodinger problems, ground state solutions, least energy nodal solu- 
tions, (nodal) Nehari set, mountain pass algorithm. 

The author is partially supported by a grant from the National Bank of Belgium and by the program "Qual- 
itative study of solutions of variational elliptic partial differential equations. Symmetries, bifurcations, singular- 
ities, multiplicity and numerics" of the FNRS, project 2.4.550. 10. F of the Ponds de la Recherche Fondamentale 
Collective. 

1 



2 



CH. GRUMIAU 



an orthogonal symmetry with respect to the hyperplane). When it is the case, we say that 
the solution respects the symmetries of V . When V = 0, this type of questions has already 
been studied. First, on the square in dimension 2, we can mention a result of G. Arioli and 
H. Koch (see [ IJ). They proved the existence of a positive symmetric ^""-function w such 
that —Am = wu^ possesses a non-symmetric positive solution (with 1 as Morse index). The 
same kind of result has also been obtained for a solution with 2 as Morse index. The proof 
is partially computer-assisted. Second, in collaboration with D. Bonheure, V. Bouchez, 
C. Troestler and J. Van Schaftingen (see [3] |7I), we proved for p close to 2 that the 
symmetries are related to the symmetries of eigenfunctions of —A. We generalize here the 
technique at some non-zero potentials V and we make numerical experiments to illustrate 
it. 

In Section [2] by denoting A,- (resp. £,) the distinct eigenval ues ( resp. eigenspaces) of 
—A + V with DBC in Hq (£2), we prove the following Theorem 1 1 . 1 1 For this, we assume 
that Xi is the unique principal eigenvalue, i.e. an eigenvalue with a related eigenspace of 
dimension 1 possessing an one-signed eigenfunction. We also require that eigenfunctions 
in E2 have a nodal line of measure 0. By using a maximum principle, these assumptions 
are satisfied at least when V e L+°°(£2) (see [6 |). 

Theorem 1.1. When V e L^^^{D.), £ L+°°(£2) and —A + V is positive definite such that 
X\ is the unique principal eigenvalue, for p close to 2, the ground state (resp. least energy 
nodal) solutions respect the symmetries of their orthogonal projections in Hq{Q,) on Ei 
(resp. E2). 

In particular, when the eigenspace has a dimension 1, the solutions respect the symme- 
tries of y. As we assumed that X\ is the unique principal eigenvalue, ground state solutions 
respect the symmetries of V. 

Depending on the structure of E2, some symmetry breaking exist for least energy nodal 



solutions (see Section 2.3 1. In fact, by a traditional bootstrap, a family of ground state (resp. 
least energy nodal) solutions {up)p^2 converges for "^-norm to functions in Ei (resp. E2). 
So, for l.e.n.s., Up does not respect the symmetries of V for p small when the projection 
is not symmetric in E2 (see Section |3] for an example). For larger p, it is depending on 
the case. In Section [273] we exhibit rectangles and V (such that the eigenfunctions in E2 
are symmetric) where l.e.n.s. do not respect symmetries of V for p large enough. So, the 
result[TT|cannot be extended to all p. 

In Section [3] as illustration, we implement the (modified) mountain pass algorithm 
(see im |9] [TOl ) to study the cases of V negative constant, piecewise constant or singular. 
We exhibit direct examples such that the solutions break down the symmetries of V. 

2. Main results 

The proofs are related to the technique defined in |3 |. This is why we just point out and 
explain the differences and we do not make all the details. The first result implies that the 
traditional Poincare's and Sobolev's inequahties are available for || ||^ := /j^|V p + y (•)^. 

Proposition 2.1. If —A + V is positive definite, y+ e L^I^{Q.) andV^ e L+"(£2), the norm 



u 



2 , 



: Jq\Vu\ +Vu and the traditional norm ||m||7,i :— Jq\Vu\ are equivalent. 



Proof. Using the Sobolev's inequalities on JqVu^ and as y e L^^^{Q.), 3C >0 such that 



|m||2< / |Vm|2+C / |Vm|2. 
Jn Jn 



Using the Poincare's inequalities and as y e L+°°(£2), 3C > and a real K such that 



|m||2 = £ / \Vu\^ + {l-e)\\uf + e / Vu^ 
Jn Jn 

>e / \Vuf- + {{l-e)C + eK) [ u^>e [ \Vu\^, 
Jn Jn Jn 



NONLINEAR SCHRODINGER PROBLEMS 



3 



where the last inequahty is obtained for £ small enough. □ 

Then, the proof of Theorem |l.l| is based on two main results. The first one shows that, 
for /7 ~ 2, a priori bounded solutions can be distinguished by their projections on 

2.1. Abstract symmetry. 

Lemma 2.2. There exists e > such that ;/ 1 1 a (jc) — A; | j < £ and u solves —Au + Vu = 
a{x)u with DBC then u = or Pe-u ^ 0. 

Proof. Similar as in Lemma 3.1 in 13], Poincare's and Sobolev's inequalities are adapted 



using Proposition 2.1 □ 



Then, we directly obtain our abstract symmetry result as in the proof of Proposition 3.2 
in ||3i. Let us remark that the result holds for any / and not just for 1 or 2 as stated 
in IS. We denote by B{Q,M) the ball in (Q.) centered at and radius M. 

Proposition 2.3. Let M > 0. For i G No, 3/5 > 2 such that, for p e (2,/5), ;/ Up,Vp G 
{u e B{0,M) : Pe-u ^ B(0, ^)} solve the boundary value problem with DBC —Au + Vu = 
A,|m|''^^m then PEjUp = ^fe, implies Up = Vp. 

These two results permit us to conclude as in Theorem 3.6 in f3\. 

Theorem 2.4. Let {Ga)aeE with E = Ei be a group acting on Hq (£1) such that, for g G Ga 
and u e Hq (Q.), 

g{E)=E, g{E^)^E^, ga = a and Sp{gu) = Sp{u) . 
For any M > I, 3p > 2 such that, for any family of solutions (Mp)p>p>2 C {m G B{0,M) : 
Pem ^ B(0, j^)} of the boundary value problem with DBC —Au + Vu — Xi\u\P^^u, Up be- 
longs to the invariant set of Gap where Up is the orthogonal projection PEUp. 



Theorem 2.4 can be used for any bounded family of solutions staying away from 0. To 
apply Theorem 2.4 at a family {up)py2 of ground state (resp. least energy nodal) solutions 
for the problem ([l_i, we study the asymptotic behavior when p ^ 2. We prove that the 
expected upper and lower bounds are fine if and only if X — Xi (resp. X2). In some sense, 
Ai (resp. X2) is the natural rescaling to work with ground state (resp. least energy nodal) 

is not a restriction. Indeed, by 
symmetries of ground state (resp. least 



solutions of problem ([TJ. Let us remark that this condition : 
homogeneity of X\u\' u in the equation ([T]), the symmetri 
energy nodal) solutions are independent of X. 



2.2. Asymptotic behavior. Let us denote {up)py2 a family of ground state (resp. least 
energy nodal) solutions for the problem ([TJ. We consider A = Ai for g.s. (resp. A„ the first 
not principal eigenvalue for l.e.n.s.) and E = E\ (resp. £„). 

Lemma 2.5. Concerning the upper bound, limsup^^^2ll"plP = linisupp^2 ( ffi-x'/p j — 
||m*|P where u^, £ E minimizes the limit functional : E ^ M. : u 1-^ J^u'^ — \ogu^. 

Proof. The proof is inspired by Lemma 4.1 in 131. First, we define Vp := + (p — 2)w 
where w £ Hq{Q.) solves the problem — Aw + Vw — A2W — A2M*log|M*| with Pew = 0. 
Then, we prove that the projection of Vp on Jf^p (resp. ^p) converges when p — > 2. 

Concerning least energy nodal solutions, in fi\, the result has been stated for n — 2. 
Here, let us remark that it works with E„ which is not specially Ej. We just need to ensure 
that Vp is sign-changing for p close to 2. □ 

Nevertheless, we need to assume n = 2 to obtain the lower bound. It is explained in the 
next Lemma. 

Lemma 2.6. Concerning the lower bound, ifn — 2 then liminfp^2||M;7|| > 0- 
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Proof. The proof is inspired by Lemma 4.4 in |i3J. Concerning l.e.n.s. (the argument is 
easier for g.s.), let e\ be a first eigenfunction in E\. By considering := f\J\ei 
Sp :— I — Sp, we show the existence of tp > such that Vp = tp{spU'lj +SpUp) belongs to 
^pHE^. 

Then, we prove that Vp stays away from zero using Poincare's and Sobolev's embed- 
dings, which concludes the proof. For this part, we need to require that Aj is the unique 
principal eigenvalue, i.e. n — 2. Otherwise, we should prove that Vp e (£1 ffi . . .E„^i)-^, 
which cannot be assumed. □ 



The two previous results imply Theorem 2.7 



Theorem 2.7. Assume that —A + V is positive definite and possesses one and only one 
principal eigenvalue (n ~2), £ L^/'^{£i) and G L+°°(n). If (Mp)p>2 is a family of 
ground state (resp. least energy nodal) solutions for equation ([T]l then 3C > such that 

ll"pll//> ^ C (^j;^ ^ 1 where i ~ I (resp. 2). If p„ — > 2 and ^ ^ ^o(^)' 

then (^j;^ ^ ''" Mp„ — > m* in Hq (£1), satisfies —Am* +Vu^, = A,m* and (m*) = inf{(o'* (u) : 
ueEi\{0},{d(^4u),u) =0}, where (^^. : E, ^R: u^ ^ j^u^ ~ u^logu^ . 
Remark 2.8. [(i)] 

(1) By a traditional bootstrap, (^x) ^ '" "Pn ^ "* ™pliss (x) ^ '"' "Pn ~^ "* ^' 
norm (see fT\). 

(2) If A < Ai (resp. A2), a family of g.s. (resp. l.e.n.s.) blows up in Hq (£2). If A > Ai 
(resp. A2), it goes to 0. So, a family of g.s. (resp. l.e.n.s.) is bounded and stays 
away from if and only if A = Ai (resp. A2). 

(3) By homogeneity of X\u\''^^u, the study of symmetries for only one value of A is 
enough to conclude symmetries for any A > 0. 



(4) By combining Theorem 2.4 and Theorem 2.7 we obtain that gound state solutions 



for p close to 2 respect the symmetries of their projection on £1 . As first eigen- 
functions are unique up to a constant, they keep symmetries of V for p close to 
2. 



(5) By combining Theorem 2.4 and Theorem 2.7 we obtain that least energy nodal 
solutions for p close to 2 respect the symmetries of their projection on E2. 

2.3. Symmetry breaking for least energy nodal solutions. For p ~ 2, previous results 
showed that the structure of least energy nodal solutions are related to the symmetries of 
M* verifying (^^(m*) = inf {S'^,{u) : u e £', \ {0}, (d(f*(M),M) — 0}. 

In ||3l (see Section 6), on the square and for V = 0, it is proved that if does not re- 
spect the symmetries of the rectangle, i.e. is not odd or even with resepct to the medians 
(which is numerically observed), then there exists a symmetry breaking on rectangles suf- 
ficiently close to the square and p sufficiently large. In our case, this property can be stated 
as follows. 

Theorem 2.9. Let us work on a square. If V is odd or even with respect to a median but 
does not respect this symmetry, then there exist some rectangles and p such that least 
energy nodal solutions Up does not respect the symmetries ofV. 

Moreover, as Up converges for ^-norm, we are able to directly construct V such that the 
least energy nodal solutions break down the symmetry of V. It will happen once is not 
symmetric. In the next section numerical experiments illustrate this interesting case. 

3. Numerical illustrations: non-zero potentials V 

In this section, we compute the (resp. modified) mountain pass algorithm to approach 
one-signed (resp. sign-changing) solutions (see lH |9] [TOl |5l). While it is not sure that 
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approximate solutions have least energy, all the other solutions that we have found numer- 
ically have a larger energy. So, we will assume that the approximations are ground state 
(resp. least energy nodal) solutions. We also give some level curves: 1 and 2 for g.s., ±1 
and ±2 for l.e.n.s. 

Numerically, we study p = 4. Let us remark that we always obtain the same kind of 
symmetry for smaller values of p. We work with p — 4 to illustrate that the result of 
Theorem 1 1 . 1 1 seems to hold at least for a non-negligeable interval. 



3.1. Negative constant potential on a square. As first example, we consider a constant 
V such that Ai > 0, i.e. V > — Ai where Ai is the first eigenvalue of —A. So, the required 
assumptions on V are clearly satisfied. Theorem ] 1 . 1 [ holds. In particular, concerning sym- 
metries, we obtain 

(1) for p close to 2, on convex domains, ground state solutions are even with respect 
to each hyperplane leaving £1 invariant; 

(2) for p close to 2, least energy nodal solutions on a rectangle are even and odd with 
respect to a median; 

(3) for p close to 2, least energy nodal solutions on radial domains are even with 
respect to — 1 orthogonal directions and odd with respect to the orthogonal one; 

(4) for p close to 2, least energy nodal solutions on a square are odd with respect to 
the barycenter 

Numerically, we consider —Am — = defined on the square £2 = (—1,1)^ in M^. 

First and second eigenvalues of —A are given by ^ and On the following graph, 
one-signed (resp. nodal) numerical solutions have the expected symmetries. Ground state 
solutions respect the symmetries of the square and least energy nodal solutions are odd 
with respect to the center 0. Moreover, the nodal line of the sign-changing solutions seems 
to be a diagonal, as for V = (see |3 1). 

• For g.s.: max(M) = 2.18, 
A(m) =2.54 

• For l.e.n.s.: min(M) = 
-4.61, max(M) = 4.61, 
<f4(M) = 33.21 

• Starting function for g.s.: 
{x-l)iy~l){x+l){y+l) 

• Starting function for 
l.e.n.s.: 

sin(;r(x + 1)) sin(2;r(3' + 
1)) 





3.2. Piecewise constant potential on a rectangle. As second example, V is piecewise 
constant on (0,2) x (0, 1). In (|6|, it is proved that there exists just one principal eigen- 
value. So, our assumptions are satisfied and Theorem |1.1| is available. For X — 1, p = 4 
and V{x,y) —V- .— Q when x < I (resp. V+ :— 10 otherwise), following graphs indicate 
that approximations are just even with respect to a direction. Ground state solutions are 
more or less "located" in x < 1 (the side minimizing energy) and respect symmetries of 
V. Least energy nodal solutions seem to be formed by g.s. on each nodal domains and are 
even with respect to a direction. 
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• For g.s.: max{u) = 5.98, 
(A(m) = 30.98 

• For l.e.n.s.: min(M) = 
—8.67, max{u) — 6.53, 
A(m) =76.23 

• Starting function for g.s.: 
{x-2){y- l)xy 

• Starting function for 
l.e.n.s.: 

sin(7r(x + l))sin(27i:(3' + 
1)) 

If y_ =0 and V+ — 35, we get the same symmetry for g.s. but l.e.n.s. are not symmetric. 
The mass is more or less "located" in the square defined by x < 1. So, we obtain a direct 
symmetry breaking. To minimize the energy, the difference in the potential is so large that 
it is better to locate the mass in one side of the rectangle. On a square and for V = 0, it is 
conjectured that l.e.n.s. is odd with respect to a diagonal. It explains the structure of the 
approximation. 

• For g.s.: max(M) = 6.19, 
A(m) = 33.14 

• For l.e.n.s.: min(M) = —9.8, 
max(M) = 9.7, S'4(u) — 
181.09 

• Starting function for g.s.: 
(x-2)(y- \)xy 

• Starting function for 
l.e.n.s.: 

sin(7r(x + l))sin(27i:(3' + 
1)) 

3.3. A singular potential on a ball. As last example, we study singular potentials. First, 
A = 1, /? = 4 and V{x.y) ~ . ] on the ball B{0, 1). Approximations show as expected 

that the ground state solutions are radial and the least energy nodal solutions are odd and 
even with respect to a diagonal. We obtain the same symmetry as for the potential V = 
(see G]). 

• For g.s.: max(M) = 4.15, <^4(m) 
29.9 

• For l.e.n.s.: min(M) = —6.36, 
max(M) = 6.36, A(m) = 76.04 

• Starting function for g.s.: 

COs(7r(.=c2 +3/2)0.5/2) 

• Starting function for l.e.n.s.: 

cos{n{x^ + y^f-^/2) cos,{2n{x^ + 

/)0-5 cos(7i:(x2+3;2)0-5 



Second, V{x^y) = , = on the ball B(0, 1). Ground state solutions seem to be 

y (.r-0.5)2+y2 

even with respect to a direction but are not radial. One can remark the work of the sin- 
gularity on the level curve 1. Least energy nodal solutions are just odd with respect to 
a direction. The mass is a little bit attracted by the side x < (the side minimizing the 
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energy). 




• For g.s.: ma\{u) = 4.41, S'4{u) = 
18.74 



• For l.e.n.s.: min(M) = —6.25, 
max{u) = 6.25, S'4{u) = 76.23 

• Starting function for g.s.: 





cos(7i;(x2 +3/2)0-5/2) 



y 



X 



(0.5,0) 



• Starting function for l.e.n.s.: 

cos(7i:(x2 + /)0-5/2) cos(2n;(x2 + 
3;2)0-5 cos(7r(jc2+/)0-5 
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